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Abstract 


Pressure data measured by the British Royal 
Aircraft Establishment for the AGARD SMP 
tailplane are compared with results calculated 
using the transonic small perturbation code 
XTRAN3S. A brief description of the analysis Is 
given and a recently-developed finite difference 
grid Is described. Results are presented for 
five steady and nine harmonically oscillating 
cases near zero angle of attack and for a range 
of subsonic and transonic Mach numbers. 

Nomenclature 

c airfoil chord, m 

Cp pressure coefficient 

•k 

Cp critical pressure coefficient 

C normalized unsteady pressure 

" coefficient; first harmonic of Cp 

divided by oscillation amplitude in 
radians 

c r wing reference chord, m 

f oscillation frequency, Hz 

“ c r 

k reduced frequency, •jy— 

LE leading edge value 

M freestream Mach number 

OB value at far spanwise boundary of 

finite difference grid 

r function defining Instantaneous 

position of wing surface, z = r(x,y,t) 

t time, s 

TE trailing edge value 

TIP wing tip value 

V free stream velocity, m/s 

x, y, z coordinates of a right hand Cartesian 
system with origin at wing root 
leading edge, positive x in downstream 
direction, y in spanwise direction, 
and z up, m 

a wing angle of attack, degrees 

Y ratio of specific heats 
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n fraction of semi span 

5, n, c transformed coordinates in x, y, and z 
directions respectively 

$ perturbation velocity potential 

t transformed time, t * t 

us angular frequency, 2nf, rad/s 

Introduction 

The transonic speed range is a critical 
region for many aeroelastic phenomena such as 
flutter and divergence. In the past, analytical 
methods have been unable to accurately model the 
nonlinear transonic aerodynamics and analysts 
have used linear theory for estimates of 
transonic aeroelastic behavior. Reliance was 
placed primarily on tests of scaled aeroelastic 
models and on flight tests of the aircraft for 
flutter clearance and aeroelastic deformation 
effects. Recently, considerable progress has 
been made in calculating steady transonic flows 
about aircraft using finite difference methods 
to obtain numerical solutions of the flow equa- 
tions. Significant progress is also being made 
toward developing finite difference methods for 
unsteady flows which may eventually lead to 
accurate transonic aeroelastic analyses. 

For two-dimensional flows, methods based on 
the transonic small perturbation (TSP) equation 
have been developed and extensively applied (see 
ref. 1-2 for example). These methods have been 
extended to include viscous effects 5 * 4 , non- 
isentropic effects, 5 and wing-canard 
configurations. 5 For three-dimensional flows, 
the XTRAN3S program has been developed' by the 
Boeing Company under USAF contract. It treats 
an isolated planar wing including aeroelastic 
deformation effects and unsteady motion. The 
XTRAN3S program has been implemented on the 
Control Data Corporation VPS32 computer at the 
NASA Langley Research Center. A variety of 
applications are being made in order to evaluate 
its applicability to several types of wings. 

Several organizations have expended 
significant effort to measure static and oscil- 
latory pressures on wings at transonic speeds 
for use In evaluating the computational methods 
and to Improve the understanding of unsteady 
transonic flows. For example, configurations 
tested at the Langley Research Center Include a 
clipped delta wing with a 6 percent thick 
circular arc airfoil section, 5 an advanced tran- 
sport wing with several oscillating controls, s 
an oscillating rectangular wing with a super- 
critical airfoil, 10 and a flexible supercritical 
wing from the DAST ARW-2 vehicle. 11 Pressures 
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for several of these have been calculated using 
XTRAN3S. 12 * 18 Calculations have also been 
made 14 * 1 * for the F-5 and LANN wings lb which 
were tested at NLR. 18 ' 17 The British Royal 
Aircraft Establishment (RAE) has recently 
measured oscillating and static pressures on a 
model called the AGARD SMP tall plane model. 18 
This configuration is sufficiently different 
from the others that it supplements them for 
verifying transonic codes. Comparison of some 
of these experimental results with calculations 
made using the XTRAN3S code Is the subject of 
this report. 

In this paper a brief description of the 
tallplane model and the test program Is 
presented first. Then the XTRAN3S program and a 
recently improved finite difference grid and 
coordinate transformation are described. 
Comparisons of calculated and measured steady 
and unsteady pressures for several subsonic and 
transonic Mach numbers are presented and 
discussed. 


Description of RAE AGARD Tailplane and Test 

The planform of the RAE tallplane model L/ 
Is shown in figure 1. The model was molded from 
graphite fiber and has a full span aspect ratio 
of 2.41, a taper ratio of 0.27, and a leading 
edge sweep angle of 50.2 degrees. The 
root chord, c r , is 0.572 m (1.88 ft) and the 
span is 0.442 m (1.45 ft or 0.773 c r ). The 
airfoil is approximately an NACA 64A010.2, and 
measured ordinates for the model were used In 
the present calculations. A hydraulic drive 
mechanism was used to oscillate the model In 
pitch about an axis at 68.2 per cent of the root 
chord with an amplitude of approximately 0.4°. 
The oscillation frequencies were 3, 12, 33, and 
70 Hertz, and the model was oscillated about 
steady angles of attack of up to ± 5 degrees. 
Only the dynamic data for angles of attack near 
0° and for frequencies of 33 and 70 Hz are used 
herein. 

Pressures were measured along five span 
stations (figure 1) with 20 pressure tranducers 
located along the upper surface of the wing at 
each span station. The same pressure gages were 
used both for the steady and unsteady 
measurements. The data acquisition system 
restricted simultaneous measurements to only two 
chords. The data give the static pressures, and 
the mean, the magnitude, and the phase of the 
oscillating pressures. The tests were conducted 
in the RAE 3-foot wind tunnel at RAE Bedford and 
were run for Mach numbers ranging from 0.65 to 
1.20 in the slotted transonic test section, and 
for Mach numbers 1.32 to 1.72 in the closed 
supersonic section. Tests were performed at 
constant pressure, and the Reynolds number at a 
Mach number of 0.86 was 3 million based on the 
wing semispan. Transition strips were also 
placed on the model at 0.075 chord. A more 
complete description of the model and tests is 
given in reference 18. 


, XTRAN3S' is 

M 2 (* t + 2 * x ) t = tU-H 2 )*, + F* 2 + G *yJx 

'y Vy'y v z'z (l) 

where the spatial coordinates x, y, and z are 
normalized by c r , the reference chord, and t 
Is normalized by c r /V. Here the time scale 
factor IT of ref. 7 is 1.0. The perturbation 
velocity potential $ is normalized by c r V. 

The coefficients for the nonlinear terms in 
equation (1) can be defined as either 18 

F = - 2 (y+1)M 2 

G - \ (y-3)M 2 (2) 

H ■ -(y-l)M 2 

or 2u 

F = - \ [3-(2-y)M 2 jM 2 (3) 


2 

H = - M‘ 


For F * G = H * 0 the linearized unsteady 
potential equation is obtained. The 

coefficients given by equation 2 are used for 
the calculations of this report for the RAE 
tallplane. 


The boundary conditions Imposed on the 
outer edges of the computational region in the 
original XTRAN3S were 


4 > - 0 

4>x + 4>t = 0 

4 Z = 0 

4>y = 0 

4y = 0 

Uzi = t*x + 4tJ = 0 


upstream 
downstream 
above and below 
wing root 
far spanwise 

wake 


where [J indicates jump in quantity across the 
wake here. 


The wing flow tangency condition is 

*Z = R x + R t ; z = 0± ’ X LE — x — X TE 
where R = r/c r . 

The outer boundary conditions have recently 
been Improved by Whitlow 21 by Implementing 
characteristic or "nonreflectlng" boundary 
conditions. These revised conditions were 
Incorporated into the present version of the 
program used for the calculations for the RAE 
tallplane. The revised conditions are 


XTRAN3S Program Description 


j (^jj- + 7 ^-) 4> t - 4>x * 0 upstream 


The modified unsteady small perturbation 
(TSP) potential equation that is solved by 
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7 { ' IT + T7 5 *t + *x 
7 *t + *z * 0 

* t - ♦z * 0 
*= 0 


2 *t + *y = 0 


0 downstream 
above 

below 
wing root 

far spanwise 


where 

B = 1 - M 2 + 2F* X 
0 = 2M/1 + 1/B 


Calculations have indicated that these 
characteristic boundary conditions significantly 
reduce the reflection of disturbances from the 
computational boundaries. Z1 


Coordinate Transformation 


The finite difference grid contains 

60 x 20 x 40 points In the x, y, and z 
directions for a total of 48,000 points. In 
physical space it conforms to the wing planform 
and the computational region is mapped to a 
rectangular domain using the shearing 

transformation 


low taper. The grid has been revised such that 
In the physical domain, the computational region 
Is a rectangular box thus alleviating the 
skewness of the grid. 

Versions of this type of grid arrangement 
have been given previously. 14 » ** In ref. 15 a 
smooth stretching was used to map the regions 
from the wing to the upstream or downstream 
boundaries. Good results were obtained for the 
F-5 wing which is highly swept and highly 
tapered. In ref. 14, equation (5) is retained 
on the wing and similar equations are used in 
the upstream and downstream regions with c(y) 
replaced by the local distance from the wing to 
the boundary of the region. Such a technique 
permits analytical evaluation of £ x and £y, 
but results in discontinuous values of these 
quantities at the leading and trailing edges. 
Good results have been also obtained for the F-5 
wing by that method. 14 

Herein a finite difference grid is used 
that is similar to that of ref. 15. In the 
physical domain the grid is described 
analytically to give a smoothly- varying mesh 
spacing. The values of C x and 5y are then 
calculated numerically using second order finite 
difference formulae at each point in the 

computational domain. The ^-distribution of 
points (for all values of n) in the 4-n domain 
is chosen to be the same as the x-distribution 
of points along the root chord in the physical 
domain. 


£ * t(x,y), n = y, C = z, and t * t 


Finite Difference Grid 


In computational space equation (1) becomes 
M 2 ^ [|- * t + 2* 5 J - U(1-M 2 )5 x * 5 

+ FE x Z <f> £ 2 + + 2G£y*£* n + Gf n 2 J 

C (4 

+ + + VcVt + *n )J1 

+ fe ^ + V + + VJ 


The grid used for the computations 
presented in this report is shown in figure 2 
along with the wing planform. The grid is 
defined in the following manner. First, 
equation (5) is used on the wing planform. The 
x-spacing of points along the, chord essentially 
is that developed earlier 12 . 1 ^ with 39 points 
along the local chord; 38 points are equispaced 
from x/c - 1/38 to 1.0, and an additional point 
is located at x/c = 1/(3 x 38). The grid 
extends from 20 wing root chords upstream of the 
wing root leading edge to 20 chords downstream 
of the wing root trailing edge. Eleven grid 

lines are used ahead of the wing and ten aft. 
The wing leading edge is centered between grid 
lines. The spacing ahead of the wing is given 
by 


The original version of XTRAN3S used 


£(x,y) 


x - x LE (y) 

= “ETyT 


(5) 


to transform the streanwise coordinate. From 
equation (5), S x * l/c(y), and ty can be 

calculated analytically from equation (5) and 
the equations describing the wing planform. 
Using equation (5) results in a physical region 
whose streamwise extent Is proportional to the 
local chord at each span station and for highly 
swept and tapered wings results in a highly 
skewed mesh in the far field of the physical 
domain. The resulting skewness of the grid led 
to numerical Instabilities that restricted the 
application of XTRAN3S to wings of low sweep and 


x 1 u * X LE * A 1 ’u + Vu ’ ’u " 1 ”” 1 UP 
. X UP ' X LE X 1 * X LE 

A - x - —X 

’up^up ' ^ *up " 1 

A 1 = X 1 ' X LE ' A 3 

where 

x l “ xi (y) = x-coordinate of first point 
the on wing 

X LE = X LE^ = x-coordinate of leading edge 

x„p = x-coordinate of the upstream boundary 
(constant) 
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and i u Is the Index of the point starting at 
the leading edge and running to the upstream 
boundary. Similarly downstream, 

X 1 . = X TE + B l'd + Vd* i d = l ’ ”• ’ON 

a 

n X DN X TE \ - X TE 


and 


3 i M 2 . i) 

’DN l DN 11 

B 1 = X W ‘ X TE ' B 3 


<DN * 1 


where 

x u = x (y) = x-coordinate of first 

w point downstream of wing 

x T)r = x T(r (y) = x-coordinate of the 
trailing edge 

X DN ■ x-coordinate of the downstream 
boundary (constant) 

and id is the index of the points starting 
aft of the trailing edge and running to to the 
downstream boundary. 

Fourteen rows of points are used along the 
wing (including the points inboard of and on the 
plane of symmetry). The rows are distributed 
along the span using a cosine distribution 


n j ‘ n TIP C0S ^2 ^ jw - 1^’ j = l ’ j 


jw 


where n TIp = span/c r and j w = number of 
grid rows p on the wing (14 here). The first 
spanwise station is at 12% span and the last one 
at 99.3% span. This distribution of grid points 
is used to emphasize the definition of loads in 
the tip region which are critical for 
aeroelastic analysis. 

Outboard of the tip, the wing planform is 
extended smoothly to the far spanwise boundary 
at 1.5 spans. The midchord line of the wing is 
extended to the far spanwise boundary as a 
parabola that has the slope of the wing midchord 
at the tip and is perpendicular to the far 
boundary. The wing leading and trailing edges 
are extended using cubic equations that match 
leading and trailing edge slopes and intersect 
the far spanwise boundary perpendicularly at one 
half the tip chord fore and aft of the midchord 
extension. The equations defining these grid 
points are 


x LE ( n ) = x LE TIp + V n "'TIp) + D 2^ n "’TIP^ 


where 


+ D 3 ( n - n TIp) 
dx, c 

D 1 -TT (tip 

°2 ~ 

An n 


An r 


An n 


= n 0B _ n TIP 


o T ip ■ n-coordinate of wing tip 

n nR = n-coordinate of the far spanwise 
extent of the grid 

x. F = x-coordinate of wing tip leading 
Lt TIP edge 

x Le * x-coordinate of leading edge of 
Lt 0B outboard extension of the wing 

Similarly for the trailing edge extension 
x TE ( n ) = X TE TIp + E l (n ' n TIP^ + E 2^ n_n TIP^ 2 

+ ^(n-nyjp) 

where 


dx 


"I 


TE 


dn 


TIP 


E 2 " / 2 C_2E l An 0 + 3(x TE nR 
An Q OB 


X TE ^ 
IE TIP 


E 3 - 3 “ ^( x Tr ” X TF 

J Ang- 3 1 v lb 0B ie T1P 


and 


x TF = x-coordinate of wing trailing 
l TIP edge at the tip 

x-j-p = x-coordinate of trailing edge of 
^OB outboard extension of the wing 

The outboard extension of the midchord line 
intercepts the right hand boundary at 

dx„ 

\b = Xm tip + l TtpJin o 

— SW ■ ( “> * E > ,/2 

and offsetting by ±c TIp /2, 

X LE = ” C TIP^ 2 

OB n 0B 

X TE “ + C TIP^ 2 

lb 0B M 0B llP 

The same streamwise distribution of points that 
is used on the wing is used in the outboard 


l -2D, Aru + 3(x. r ‘ X,r )J 
i u L - 0 B ll tip 

region. 


The 

[DjAru - 2(x, r x i f 

1 U Lt 0B TIP 

outboard 

equation, 
the wi ng 
located 
applyi ng 
gives 
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jf-o = number of spanwlse points outboard of 
UB wing tip 

j = Index of the points outboard of wing tip 

n 0l = n- ^ oca t1° n of the first point outboard 
of the wing tip 

In previous versions of XTRAN3S, points 
Inboard of the wing root were located by 
extrapolating the first outboard points through 
the wing root using straight lines. For the 
calculations for this wing the x-coordlnates of 
the points located inboard of the plane of 
symmetry are reflected as mirror images of the 
points at the first station outboard of the 
plane of symmetry. The values of 5 X and 5y 
are then calculated by the same finite 
difference equations as previously Indicated. 
This grid variation, first Implemented by Dr. 
John T. Batina of Langley Research Center, has 
been shown to improve the pressures calculated 
along the centerline. Further improvement has 
been obtained by Batina by slightly rounding the 
wing apex. This change was not implemented In 
the version of XTRAN3S used here. 


The grid extends 25 root chords above and 
below the wing. Twenty rows of points are used 
above and 20 rows are used symmetrically below 
the wing. The distribution of points Is that 
previously developed and applied. » 21 

The resulting grid Is smooth and embeds the 
planform smoothly. This grid appears to be 
reasonable, but since to date only limited 
variations have been investigated for three 
dimensional configurations, the grid cannot be 
considered an optimal choice. The grid Is 
reasonably fine on the wing but Is relatively 
coarse off the wing. 


Results and Discussion 


RAE AGARD Tail plane Model 


Calculations have been made with XTRAN3S 
for the RAE tailplane near zero angle of attack 
and for Mach numbers of 0.65, 0.80, 0.86, 0.90, 
and 0.95. The numerical algorithm used In 

XTRAN3S Is an alternating-direction Implicit 
scheme but with some terms treated explicitly. 
Treating some terms explicitly leads to a CFL 
(Courant-Frledrlchs-Lewy)-type stability 

restriction In the time accurate finite 

difference solution. For the calculations for 
the tailplane model, this stability limit was 
found to be near at = 0.0075. This small time 
step results In a large number of steps per 


number of steps per cycle Increases with 
decreasing frequency. Calculations were made 
only for freqencles of 33 and 70 Hertz which 
required approximately 2000 and 1000 time steps 
per cycle respectively. Generally three cycles 
of oscillation were calculated starting from a 
converged steady flow solution. The last cycle 
was analyzed for the Fourier components. For 
the steady calculations the program was run In a 
time asymptotic manner for 1500 time steps to 
assure convergence. On the VPS32, XTRAN3S 
requires approximately 1.5 seconds of CPU time 
per time step. The code Is not vectorized 
except for that generated by the compiler. 

The measured airfoil ordinates have been 
used to define the wing surface slopes required 
for Input to XTRAN3S. The measured ordinates 
were fitted with a parametric spline with 
smoothing, 23 and the slopes were calculated from 
the spline fit. 

The experimental data points presented 
herein were obtained from a digital tape 
supplied by the RAE and are presented unedited 
In the figures. 

Results for Steady Pressures .- The variation 
of the calculated static pressure with Mach num- 
ber is shown in figure 3. The pressures for the 
upper surface (-Cp) are shown for an angle of 
attack of - 0.3 degrees. Typical of highly swept 
wings, the leading edge pressures change from 
compression at the root to expansion at the 
tip. The pressures are supercritical at 

M = 0.86, and at M = 0.90 a strong shock wave Is 
formed which Is particularly pronounced near the 
tip. At M = 0.95, the strong shock has moved 
rearward to near the trailing edge. 

The calculated and measured pressures are 
compared in figures 4a-4e. Points for repeated 
measurements are shown with plus signs within 
the symbols. The experimental data for this 
wing indicated that zero bending moment was 
obtained near an angle of attack of 0.3 to 0.4 
degrees which resulted from tunnel flow 
angularity. The results herein (both steady and 
unsteady) were run at -0.3 degrees angle of 
attack, which does not correspond to the experi- 
mental estimate. However, some indication of 
the effects of a small angle of attack can be 
obtained from the calculated results. The cal- 
culated lower surface pressures are more nearly 
comparable to the experimental upper surface 
pressures as a result of this sign difference. 

For M = 0.65 the agreement is reasonable 
for the inboard stations (fig. 4a). However, 
the agreement deteriorates at the outboard sta- 
tions. The pressure expansion over the forward 
portion of the airfoil is generally underpre- 
dicted. Similar trends are noted also for M = 
0.80 and M = 0.86 but with Improved agreement 
Inboard (figs. 4b and 4c). There Is a strong 
peak In the measured pressures near the tip 
leading edge for these two Mach numbers which 
may indicate a viscous effect around the sharply 
cut-off tip. For M = 0.90 (fig. 4d) there Is a 
shock which is particularly strong near the tip 
and whose strength is overpredicted by XTRAN3S. 
Again the agreement on the forward portion of 
the airfoil is resonable and underpredicted by 
XTRAN3S. At M = 0.95 (fig 4e), the shock 
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location Is predicted to be significantly 
rearward of the experimental location. 

Overall, XTRAN3S predicts the general 
trends In the static pressures fairly well, but 
with poor agreement near the tip. Best agree- 
ment Is obtained near critical free stream 
Mach number. An underprediction of the 
expansion pressures Is Indicated with an over 
prediction of the strength of transonic shock 
waves. One limitation of TSP theory as 
Implemented In XTRAN3S Is that only the 
streamwlse velocities are considered In the 
expressions for pressure. For highly swept 
wings, spanwlse flow effects may need to be 
considered especially for nonzero angles of 
attack. 

Unsteady Pressures .- The comparison of 
calculated and measured unsteady pressures for 
the 70 Hz cases for each of the Mach numbers Is 
shown In figures 5a-5e. The data are presented 
along each chord In the form of real and 
Imaginary components of the first harmonic of 
the unsteady pressure. For M = 0.65 and 0.80 
the calculated results are shown only for the 
upper surface; the results for the lower surface 
were Identical. For the higher Mach numbers, 
the lower surface pressures are shown with their 
sign changed. Only small differences are shown 
in the calculated upper and lower surface 
pressures for M = 0.86, with modest differences 
shown In the region of the shock at M = 0.90 and 
M = 0.95. For M = 0.65, 0.80, and 0.86 there Is 
good agreement between calculated and measured 
trends, and good quantitative agreement of the 
pressures over the inboard sections where the 
agreement for the pressures results was good. 
As might be anticipated, the overall agreement 
at the tip Is not good where the steady results 
were not very good. In addition, the XTRAN3S 
results do not show the sharp pressure peak near 
the leading edge that is shown by the 

experimental results. A calculation for the M ■ 
0.86 case has been made using the Euler 
equations (ref. 24). Although no detailed 
comparisons have been made, the trends are 

comparable, but with somewhat better agreement 
near the leading edge. 

For M = 0.90 and 0.95 and 70 Hz, the 

agreement In trends and magnitudes on the 

forward portion of the airfoil Is comparable to 
that at the lower Mach numbers. However, near 
the shock, the calculated results deviate 
significantly from the measured values. The 
unsteady results are seriously affected by the 
differences In steady flow shock location and 
strength. 

Where repeat points are shown, the 

agreement of the two measurements is very good 
(figs. 5a and 5c for example) and the symbols 
with and without plus signs almost overlay. 

The comparison of calculated and measured 
results for the 33 Hz cases Is shown In figures 
6a-6d. The M = 0.65 case was not run for this 
frequency. For 33 Hz, the value of k varies 
from 0.234 at M = 0.80 to 0.201 at M = 0.95. 
Comparison of the 33 and 70 Hz data shows about 
the same trends and magnitudes with the 
imaginary part for 33 Hz being about one-half of 
that for 70 Hz as expected. The level of 


agreement Is also comparable to the higher 
frequency results but with better agreement for 
the Imaginary parts. For M = 0.90 and 0.95 the 
calculated peaks near the shock are larger than 
for the higher frequency cases. Such peaks are 
not observed on the experimental data. 

Several cases have been calculated for the 
F-5 planform using XTRAN3S. ll *» 1:> The agreement 
with the experimental data Is better than for 
the RAE tallplane. However, the calculated 
shock wave strength Is also overpredicted for 
the F-5 for M = 0.95. The F-5 wing Is 
significantly thinner as It has a slightly 
cambered 4.8 percent thick section and should be 
more consistent with TSP theory. Applications 
of XTRAN3S to the calculation, of pressures on 
thick supercritical wings have also been 
given. lu , r2 , 14 , l/ Generally the correla- 
tion of calculated and experimental results is 
less satisfactory than for the F-5 or the tail- 
plane configurations. Improvements In the 
theory or further development of the code may 
overcome some of these current limitations. It 
should also be noted that significant questions 
have been raised concerning the accuracy of 
isentroplc potential theory 211 for two dimen- 
sional flows. A nonisentropic small disturbance 
theory has recently been developed to alleviate 
some of the problems In two dimensional codes. s 


Concluding Remarks 

An overview of the transonic small 
perturbation code XTRAN3S and a recently 

developed finite difference grid have been 
given. Data measured by the British RAE for the 
AGARD SMP tallplane have been compared with 
results calculated using XTRAN3S. Several 
steady and higher frequency unsteady cases with 
near zero angle of attack and for subsonic and 
transonic Mach numbers are treated. 

Good overall trends were obtained with 
XTRAN3S for the RAE tallplane model but poor 
agreement was found near the tip, the leading 
edge, and for strong shock waves. Extensive 
computer resources are required to run XTRAN3S 
for an Investigation of this type. Further 
effort is needed to reduce the resources 
required by increasing the numerical stability 
to permit large time steps. Further development 
is needed to treat thick supercritical wings 
with shock waves. 
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Fig. 2 Finite difference grid in the plane of 

the wing - near field and far field views. 
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(e) M = 0.95. 
Fig. 4 Concluded. 



Fig 5 Measured and calculated unsteady pressure distributions for tallplane model 
oscillating at approximately 70 Hz. 
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(b) M ■ 0.86, k « 0.222. 


Fig 6 Measured and calculated unsteady pressure distributions for tallplane model 
oscillating at approximately 33 Hz. 
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